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ABSTRACT 

Understanding the origin of AGN absorption line profiles and their diversity could help to explain the physical 
structure of the accretion flow, and also to assess the impact of accretion on the evolution of the AGN host 
galaxies. Here we present our first attempt to systematically address the issue of the origin of the complexities 
observed in absorption profiles. Using a simple method, we compute absorption line profiles against a contin- 
uum point source for several simulations of accretion disk winds. We investigate the geometrical, ionization, 
and dynamical effects on the absorption line shapes. We find that significant complexity and diversity of the 
absorption line profile shapes can be produced by the non-monotonic distribution of the wind velocity, density, 
and ionization state. Non-monotonic distributions of such quantities are present even in steady-state, smooth 
disk winds, and naturally lead to the formation of multiple and detached absorption troughs. These results 
demonstrate that the part of a wind where an absorption line is formed is not representative of the entire wind. 
Thus, the information contained in the absorption line is incomplete if not even insufficient to well estimate 
gross properties of the wind such as the total mass and energy fluxes. In addition, the highly dynamical nature 
of certain portions of disk winds can have important effects on the estimates of the wind properties. For ex- 
ample, the mass outflow rates can be off' up to two orders of magnitude with respect to estimates based on a 
spherically symmetric, homogeneous, constant velocity wind. 

Keywords: Accretion, accretion disks — Black hole physics — Line: formation — Line:profiles 



1. INTRODUCTION 

Accretion disk winds are among the most promising physi- 
cal mechanisms able to link the small and the large scale phe- 
nomena in active galactic nuclei (AGN) and to shed light on 
the physics of mass accretion/ejection around supermassive 
black holes (SMBHs). Such winds are currently directly ob- 
served, as blueshifted and broadened absorption lines in the 
UV and X-ra y spectra of a substantial fraction of AGN (e.g., 
ICrenshaw eTal.l i2003). 

In the UV band we observe, with decreasing width of the 
absorption troughs, broad absorption lin es (BALs) in about 
15 % of optically selected AGN (e.g., iKnigge etall I2008t 
iGibson et all200l : mini-broad absorption lines (mini-BALs) 
and narrow absorption lines (NALs) in about 30 % of AGN 
dOanguly & Brotherton 2008). Once corrected for selection 
effects, the intrinsic fraction of AGN hosting B ALs is es- 
timated to be ~ 30 - 40 % (lAllen et all 1201 lb : summing 
also the corrected fraction for mini-BALs and NALs, the in- 
trinsic fraction of AGN hosting blueshifted UV absorption 
lines results to be as high as ~ 70 % (Hamann et al. 2012). 
In all these cases, the absorbers are associated to resonant 
transitions of moderately ionized elements like Mg n, Al iii, 
C IV, Si v, and can be outflowing with velocities as high as 
~ 0.2 c, the main difference among them being in the width of 
their absorption troughs. Given their large width, BAL struc- 
tures are the easiest to be identified even in low-resolution 
UV spectra, and have therefore been studied in much more 
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detail than mini-BAL and NAL structures. The shape of 
BAL troughs is quite complex, showing a variety o f pro- 
files among diff^erent sources (e.g.,[T'u i'nshek 1984; Folt z et al.l 
\m% ITurnsheketai] [T988t IHall et al.i ^002). When high 
spectral resolution, high signal-to-noise observations have be- 
come available, generally also mini-BALs and NALs showed 
a variety of profiles ( e.g.,'Misawa et al.ll 2007l: IWu et al.ll2010l: 
iHamannetaTI 120111: Rodriguez Hidalg o et al.1 1201 Ih . Re- 
cently a number of AGN have been shown to display strong 
variability in their UV absorption troughs, including emer- 
gence of BAL structures and transitions from mini-BALs 
to BALs (e.g.. Ma 2002; Hamann etal.ll2008t iLeighJxeLal] 
120091: iKrongoid e t al. 2010; Vivek et al. 2012). 

In the X-ray spectra of ~ 50 % of AGN, we observe also 
lower- velocity (from ~ 100 to ~ 1000 km s"') outflow- 
ing ionized gas (the so-called warm absorber) in the transi- 
tions of joni zed elem ents s uch as C vi, O vii, O viii. Si vi 
(e.g., iMcKernan et al.ll2007l) . Higher velocity absorbing gas 
in the transitions of Fe xxv and Fe xxvi, with blueshift in 
the r ange 0.003 - 0.3 c i s observed in a number of AGN 
(e.g., 'Poun ds et al.l 1200?. 'Risaliti et al.' '200 l ICappil l2006t 
Miniutti et al.l l2007l: TTurner et al. 2008; Giust ini et aLlboi lb ; 
their fraction among the AGN population is currently un- 
known, the only reliable estimate bein g of about 30 - 40 % 
for a sample of low redshift AGN (iTombesi et al.l j^OlO). 
X-ray BALs outflowing up to 0.7 c ha ve been also ob- 
served in a handful o f AGN (IChartas et al.l l2002. 2003, 2009] 
iLanzuisi et al.ll20T2b . The ionization state of the X-ray ab- 
sorbers spans a large range of values, from e.g. C vi in 
the case of warm absorbers up to Fe xxvi for X-ray BALs. 
The observed high-velocity X-ray absorbers profiles are quite 
complex and display strong temporal variability both in depth 
and velocity shift; also for the case of the lower-velocity X- 
ray warm absorbers, deep observations have revealed sig- 
nificant complexities in the absorption trough profiles (e.g.. 
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Figure 1. The three "typical" streamhnes for accretion disk winds. 1: MHD- 
typical, polar and convex streamline. 2: hydrodynamical-typical, equatorial 
and concave streamline. 3: a radial streamline typical of spherically symmet- 
lic winds. 

iRamfrez et ani2005l: [Pounds & VaughaDll201 Ih and variabil- 
ity in ionization and v elocity shift (e.g., see iLonginotti et al] 
l2010l:lMatt etal.ll201lL for some recent results). 

Velocity and ionization state are only two of the physical 
properties of AGN winds; the mass and energy fluxes car- 
ried by them are other fundamental properties that we wish to 
be able to measure in order to quantify the actual impact of 
these winds on the surrounding environment, i.e., the amount 
of feedback. To measure the mass and energy flux, the den- 
sity of the wind material and the distance between the cen- 
tral SMBH and the place where absorption occurs must be 
known. Unfortunately, these two quantities are difficult to 
measure, and estimates for them have always relied on several 
assumptions. In particular, a geometrically and optically thin, 
constant-density, spherically symmetric, single-zone of gas in 
ionization equilibrium and outflowing with a uniform velocity 
has been usually assumed to derive constraints on the distance 
of the absorbing material from the central ionizing source. 

Using these assumptions, reported estimates of the dis- 
tances of absorbers from the central SMBH range from the 
inner regions of the accretion disk for the high-velocity X- 
ray BALs (e.g., Chartas et al. 2009), to the parsec -s cale torus 
for the X-ray warm absorbers (e.g., IBlustin et alJ I20d5h. to 
the k iloparsec-scale for some UV BALs ("e.g.. IPunn et alJ 
120101) . Currently the uncertainties on the distances are rel- 
atively large, which translates to a high uncertainty on the 
mass outflow rate and on the kinetic energy injection associ- 
ated with such winds (e.g. JReeves et alJ2003l : lKrongold et alJ 
12007; Cappi et al. 200%. 

Theoretical arguments and numerical simulations have 
been used to show that accretion disks are able to launch and 
accelerate powerful outflows by various physical mechanisms 
such as radiation pressure, magnetic pressure, and thermal 
pressure (see e.g. Konigl 2006; Proga 2007; Everett 2007, for 
reviews on the subject). The driving mechanism and the ac- 
creting system physical parameters (e.g., accretion rate, black 
hole mass, magnetic field configuration, UV/X-ray flux ratio) 
determine the wind characteristics (e.g.. iPro^ 120051) . De- 
spite the many parameters that combine themselves in com- 
plicated ways, only a few geometries of streamlines are shown 
by AGN accretion disk wind models. 

Figure [T] shows the main general shapes of streamlines that 
are predicted by models of AGN accretion disk winds. Num- 
ber 1 is a convex, po lar streamline, typical of magnetic driv- 
ing scenarios (e.g., Blandford & Payne 1982; Lovelace et al.' 
[T99h.Krasnopolskv et al.,2003; Anderson et aL2005; .Everett. 



l2005h . Number 2 is a concave, more equatorial 
streamline, typical of Une-driven (e.g., iProga et al] 120001 ; 
iProga & Kallman' '2004*) and thermally dri ven disk w inds 
(e.g'. Woods et al. 1996; Font et al. 2004; LuketiceLalBoiO). 
What is usually assumed in order to deduce quantities such as 
mass outflow rate, ionization state, or absorber distance from 
the continuum source is the configuration number 3, i.e. a 
radial flow. 

However, the different geometries correspond to very dif- 
ferent wind properties: the mass flux m„„f across the area of a 
flow tube, A(r) is 

rhout ^ A{r)u(r)p(r) (1) 

where v is the velocity normal to the area, and p is the density. 
For a simple radial wind with a constant velocity, A(r) oc 
and p(r) oc r"^, whereas for a cylindrical wind (with constant 
velocity) A(r) = const and p(r) = const. This density depen- 
dence on the flow geometry has received little consideration in 
previous works on estimating th e flow mass and energy rates 
using observations (however, see IWaters & Proiall2012[) . 

In this article we will report on results from our first at- 
tempt of systematically addressing the following questions: 
what drives the appearance of different line profiles? Can 
we distinguish among different geometries for accretion disk 
winds from observations? How reliable are estimates of wind 
properties, when a radial or spherical flow is assumed? The 
secondary questions that will lead to answer the former ones 
are then: what are the implications of radial flow assumptions 
when deriving physical quantities related to the wind? Do 
different wind physical characteristics produce different line 
profiles? If yes: how and why? Our method to compute line 
profiles is introduced in Section 2; results are presented in 
Section 3 and discussed in Section 4; finally, conclusions are 
presented in Section 5. 

2. METHOD 

We consider several examples of disk wind models, with 
an increasing level of completeness or complexity or both. 
We present results in terms of synthetic absorption line pro- 
files predicted by these models. Specifically, we first con- 
sider a simplified, isothermal, thermally driven disk wind 
model and use steady state solutions as pr edicted by hydro- 
dynamical simulations similar to those of iFont et alj (l2004t) 
in Section 13.11 t hen we cons i der hy drodynamical simula- 
tions as given by iLuketic et al.l (120101) . where the effect of a 
central X-ray illuminating source is explicitly taken into ac- 
count, showing t he e ffects of varying ionization state across 
th e flow (Sectionl3.2l); fina lly we consider the solutions given 
bv lProga & KallmanI (12004 ). showing the effects of the depar- 
ture of the flow from a steady configuration (Section [33] ). 

Our method to compute absorption line profiles adopts a 
radiation transfer treatment that is simplified to the bare min- 
imum, nonetheless capturing the most important physical and 
geometrical effects, and as we will show is useful to explore 
the effects of different wind geometries. Specifically, we ig- 
nore scattered radiation and we assume that the continuum 
photons are emitted by a point source. Therefore we deal with 
proper "line of sight" (LOS), instead of "cylinder of sight" and 
line profiles observed at infinity are sensitive only to the ra- 
dial component of the wind velocity only along a given LOS. 
This kind of approach is viable for analyzing e.g. X-ray prop- 
erties of AGN winds because of the compactness of the X- 
ray source continuum. Having the distributions of the density 
p and the radial velocity Vr for each wind model, we com- 
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pute synthetic absorption line profiles corresponding to dif- 
ferent LOS or equivalently inclination angles. Our procedure 
is the following: First, we search for resonance points given 
the continuum frequency, Vg along a given LOS using the res- 
onance condition 

Vo C 

where v is the Doppler shifted flow frequency with respect to 
the continuum, and v,- is the velocity along the LOS. We then 
use the Sobolev approximation to compute the optical depth: 

T = —p ., . w , I (3) 
V \avr{r)/ar\ 

where the line opacity /c is a product of the oscillator strength 
and the abundance of a given element and ionization factor 
for a given ion, and dvr(r)/dr is the velocity gradient along 
the LOS (with our approximation it is just the gradient of the 
radial velocity). To focus on the effects of the density and ra- 
dial velocity distributions across the flow, we set kc/v = 1 and 
compute the optical depth r,, , at a given resonance point /, and 
the total optical depth, Ty at a given frequency by summing up 
the optical depths at all resonance points that we found 

Tv = ^ Ty,i (4) 
/ 

We finally solve the radiative transfer equation: 

/ = loe-''. (5) 

Many already studied how synthetic line profiles depend 
on the wind geometry and viewing angle. However, most 
of previous studies used kinematic models of accretion 
disk winds a nd were applied to cataclysmic variables (e.g. , 
Drewl[T 987^ 'Shlos man & \^teIiol [T99l iKnigge et al.l [l99l 
Long & Knigge 2002^! These models were axisymmetric. 
They assumed streamlines to be straight lines and the wind 
velocity to be some prescribed function of the position along 
the streamline. By varying one of the free parameters (namely 
the inclination angles that the streamlines made with respect 
to the disk plane), such models permit the exploration of var- 
ious wind geometries, but only to some extent. One of the 
results of these studies is that the formation of blueshifted 
broad absorption troughs (i.e., P Cygni profiles) can be nat- 
urally obtained for a variety of wind geometries. Our work 
is an extension of those earlier studies as we investigate how 
the absorption troughs depend on the geometry of the stream- 
lines, the opacity, and the ionization state but instead of using 
kinematic models we use solutions of hydrodynamical equa- 
tions obtained by numerical simulations. More sophisticated 
calculations of wind photoionization and radiative transfer 
are possible, however they are relatively time consuming and 
there fore less suitable t o the exploration of various geometries 
(e.g.. lSim etal.ll2010bl) . 

3. RESULTS 

3.1. Effects of the Overall Wind Geometry 

We first consider 2D, axisymmetric, isothermal, steady 
state disk w ind models. We use models computed by 
iLuketic et al.l (I201Q) and refer a rea der to this work for d etails 
about the simulations. We note that lLuketic et al.l (120101) per- 
formed these simulations to test their code and setup against 
results presented bv .Font et al.. (.2004) . 



The latter showed that a relatively wide range of the wind 
geometries and shapes of the streamlines can be modeled by 
varying the slope of the radial profile of the density along the 
disk atmosphere (the base of the wind). In particular. Fig- 
ure |2] shows examples of three different wind geometries as- 
suming a power-law density profile: p{r) oc r^ and three val- 
ues of a: 1. 5 (model A), 2.0 (model B), and 2.5 (model C). 
As found bv lFont et al.l (120041) . all the three models predict a 
very similar shape and divergence (i.e., A{r)) for the stream- 
lines originating from the inner most part of the disk; for other 
streamlines (or middle and high inclination angles), some dif- 
ferences and trends among various models start to appear The 
flattest radial density profile along the wind base (model A) 
corresponds to a relatively polar geometry, while the steepest 
profile (model C) is the one that most resembles the spherical 
outflow case (quasi radial, with spherical divergence). 

We probed the winds along five different LOS, indicated 
with the colored dashed lines in Fig. |2l the radial velocity, 
density, radial velocity gradient, and finally optical depth are 
plotted in Fig. |3]for each model and for each inclination an- 
gle. All the quantities are arbitrarily normalized, as at this 
point, we simply want to explore the qualitative behavior of 
the different geometries: to this end, what is relevant is the 
relative contribution of the different parts of the flow and not 
their absolute magnitude, the latter being strongly dependent 
on the actual details of the system such as the central mass, 
the Eddington ratio, and the spectral energy distribution. 

Figure |4] shows the resulting absorption line profiles using 
the same color code as in the previous figures. To show the 
profile shapes in more detail, the bottom panels of Figure |4] 
show results where the opacity has been multiplied by a fac- 
tor so that the point in velocity space where the opacity is 
maximum reaches the zero flux level for each LOS. 

The most obvious difference between the profiles predicted 
by winds of different geometries is in the shape of the blue 
wing: very sharp for model A and very smooth for model C. 
Model B is an intermediate case: similar to model A at high 
inclination angles (/ ~ 80°), to model C at middle and low 
inclination angles (/ < 60°). The main effect of / on the pro- 
files, regardless of the wind model, is on the profile width: the 
larger the inclination the narrower the profile. However, the 
profile sensitivity to / is model dependent. In particular, the 
profile width is most sensitive to / in model A and the least in 
model C. This is simply a reflection of the fact that model A is 
polar whereas model C is almost spherically symmetric. The 
right panels of Figure [3] show that for Model C, at large radii, 
the flow properties do not depend on / (profiles for various / 
collapse onto one power-law) as expected for a spherical solu- 
tion, whereas the panels in the left and middle columns show 
that the flow properties are strong functions of /. 

For the highest inclination angle (light blue lines in Figure 
[3]) in Model A and C, the radial profiles for the opacity and 
the radial velocity gradient appear to be noisy. In the case 
of Model A, this is a numerical effect which is related to the 
nearly constant velocity, that makes its gradient very small; as 
for Model C, the noise reflects the fluctuations in the radial 
velocity profile close to the base of the wind (i.e. the disk 
surface). 

In this subsection, we consider the isothermal disk wind so- 
lutions. One of the key consequences of assuming constant 
temperature is that the wind velocity along a given streamline 
is an ever increasing function of the radius. In other words, 
those models have infinite acceleration zones. This can be 
seen most clearly in the top right panel of Figure [3] showing 
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Model A 





Model C 




Figure 2. Streamlines for the three simple isothemial disk wind models adopted from lLuketic et al.lj2010h (see also lFont et al.l2004h . The models were computed 
assuming a density profile at the wind base as the following function of radius, ro'- p(r) k r^" for a = 1.5,2.0, and 2.5, from left to right. The colored dashed 
hnes represent the five LOS for which we calculate absorption line profiles: = 15, 30, 45, 60, and 80° for the black, red, green, dark blue, and light blue line, 
respectively. The length scales on the figure are in units of the radius at which the gas internal energy equals the gravitational energy of the central object located 
at (0,0). 
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Figure 3. Radial profiles for the velocity, density, radial velocity gradient, and finally line optical depth, for model A, B, and C from left to right. The color code 
is as in Figure|2] Radii ai'e in units of the radius at which the gas thermal energy equals the gravitational energy of the central object, while velocities are in units 
of sound speed. 



the radial velocity for model C and various /. In model C, 
the flow is almost radial so that the velocity plotted in Fig- 
ure[3](the radial component) is very similar to the total veloc- 
ity. Note that this is not the case for the plotted velocity of 
model A that could give an impression that the velocity satu- 
rates at a finite level. This saturation is not caused by a finite 
acceleration but by projecting the total velocity onto a LOS 
(in effect, taking only the radial component of the total veloc- 
ity). This projecting significantly affects the plotted velocity 
for model A and to some extent for model B. 



A more detailed analysis shows that for the case of polar 
flows (model A), the optical depth is a weak function of ra- 
dius at large radii. This explains a relatively flat line profiles 
at large velocities and a sharp edge corresponding to the ve- 
locity at roui- The flatness of the r profile is caused by the 
relative flatness of the density and radial velocity profiles. In 
the case of more equatorial flows (model B and C), the den- 
sity scales as expected for the spherical flow like r"^ and t is 
a decreasing function of radius at large radii. Therefore, the 
line absorption reaches a maximum at some intermediate ve- 
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Figure 4. Absorption line profiles for the for each model, along the five LOS color-coded as in Figure|2] Top panels: profiles computed using the actual opacity 
of the models; bottom panels: profiles computed increasing the opacity for each LOS until the zero flux level is reached at least in one point in the velocity space. 






locity and then it gradually decreases with increasing velocity 
with the signature of the finite size of the computational do- 
main (the sharp edge) being weak or undetectable. As we 
will illustrate in next section, sharp edges or features in the 
line profiles can be produced not only as artifacts of the finite 
computational domain but also as physical consequen ces of 
the internal structure of the wind (see also the work of lDrewl 
[1987'; 'Shlo sman & Vltelio lTQQ?: Knigg e et al.l[T995h . 

Although the profile width depends on its overall shape 
is not very sensitive to / for these simple isothermal models: 
the diff'erences between profiles for different wind geometries 
seen at different inclination angles are not very significant - 
most profiles are flat, featureless, and differ only in the posi- 
tion of the blue edge. In particular, looking at the bottom pan- 
els of Figure |4] it is really hard to distinguish among model 
B and C, and also among the different inclination angles. For 
example, the line profile for / = 60° of model B is very simi- 
lar to the profile for / - 80° of model C. Looking at model A, 
the profile shapes are now undistinguishable among different 
LOS, the only difference being in the maximum blueshift ve- 
locity and the width of the absorption trough. Given the very 
similar characteristics of the synthetic absorption line profiles 
resulting from different inclination angles and different ge- 
ometries, we conclude that different wind geometries are not 
the main contributors to the very different shapes of observed 
absorption lines. 

3.2. Effects of Finite Acceleration Zone and of Position 
Dependent Ionization State 

In the previous section, we found that as simple and easy to 
understand the isothermal models are, they are of a limited 
usage and could be a source of artifacts in the line profiles. 
For example, the position of the blue edge in the profile for 
model A would shift to higher velocities if the outer radius of 
the computational domain were increased. This shift would 
be a consequence of the fact that the computation would cap- 
ture a faster part of the wind solution. For these and other rea- 
sons, we consider next a model for which the assumption of 
isothermal flow was relaxed and instead the gas temperature 



was computed by introducing physically consistent radiative 
cooling and heating processes. We use o ne of the sim ulations 
of the thermally driven disk wind from Luketic et al.l (120100 . 
Specifically, we use a rerun of their fiducial run (C8) on a grid 
with an outer radius r - lO Orc where rc i s the C ompton ra- 
dius. For most of their runs. iLuketic et al] (1201 Oh used a grid 
with a smaller outer radius (r = 20rc). However, to examine 
and confirm the self-similarity of the flow at large radii, the 
performed several simulations on the bigger grid. Figure 
reports the geometry of the wind and the five probed LOS in 
its top panel, and the corresponding absorption line profiles in 
its bottom panels. The difference in profiles with respect to 
the previous case is striking. In this case, the inclination an- 
gle / plays a fundamental role in shaping the absorption line 
profiles. Generally, the profiles are no smooth nor broad any- 
more, and show two distinct components for a range of /. 

To associate the line profile characteristics with the wind 
properties, Fig.|6]shows the radial profiles of the radial veloc- 
ity, the density, the radial velocity gradient, the optical depth, 
and the ionization parameter for five /'s. The main differ- 
ence compared to the previous model is that the radial velocity 
reaches a maximum or two within the computational domain 
for most / and after reaching the maximum, the radial velocity 
can significantly decrease with increasing radius (see the top 
panel of Fig.|6]l. This results in large velocity gradients and in 
turn sharp maximums in r. 

However, the line forming region can be smaller than the 
size of the wind, so that only one of the line components can 
be produced by a real system. This can occur in winds where 
changes in the flow ionization limit the formation of a given 
line to a relatively small region. To illustrate this point, we 
will scale the ion fraction with the photoionization parameter 
^ - L/Annr^ where L is the ionizing luminosity, n is the num- 
ber density of the gas, and r is the distance from the contin- 
uum source. Only in special cases ^ is constant, for example 
when n oc r^^ and the flow is radial and of constant veloc- 
ity. For a radial flow with variable velocity, ^ oc v(r) so ( 
increases with increase outflow velocity (e.g.. iRamfrez et all 
120051) . However, for most /, in all the disk wind models we 
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Figure 6. Radial profiles for the velocity, density, radial velocity gradient, 
line optical depth, and ionization parameter, for the rerun of the Luketic et al] 
{2010) model. The color code is as in Figure [5] the quantities have been 
normalized to their maximum value among all the LOS. 
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Figure 5. Top panel: streamlines for the disk wind model of ILuketic et alj 
{2010); bottom panels: absorption line profiles for dift'erent inclination an- 
gles, wit h the same color co de as in the top panel, for the rerun of the C8 
model of lLuketic et alJi2010D . The velocity scale has been normalized to the 
maximum outflow velocity v„ax = 6.6 X 10^ cm s"'. 



discuss here, n decreases with r slower than r"^ and could 
even increase with r. Consequently, the photoionization pa- 
rameter decreases with increasing radius for a given LOS. 

In Figure]?] we show the effects of the varying photoioniza- 
tion parameter through the flow, on the absorption line profile 
for i - 60°. Specifically, we recalculated the profile assuming 
that the number density of a given species is a function of ^: 
we used two thresholds with a difference in £, of about two 
orders of magnitude (e.g., to mimic the difference between 
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Figure?. The effect of varying ionization parameter for the i = 60° LOS: 
in black the lower ionization state species, in red the higher ionization state 
species, in blue dashed line the total line profile without accounting for ion- 
ization parameter distribution effects. 
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O vm and Fe xxv). The profile without effects of ^ is plotted 
using a blue dashed line (it is the same profile as in Figure |5] 
plotted using a blue solid line, and it is double-peaked reflect- 
ing the two maxima in the opacity distribution along this LOS 
visible in Figure |6]l, while the profiles with the effects of the 
ionization state are: the red line refers to a highly ionized 
phase and the black line refers to a lowly ionized one. After 
taking into account the ionization state distribution across the 
flow, one can easily lose the double-peakness and find a sin- 
gle absorption trough for each ionic species. Furthermore, as 
a consequence of the radial distribution of ^ through the flow, 
one can have a greater blueshift for the lowly ionized species 
than for the highly ionized ones. This is an opposite trend to 
the one expected for a radial outflow! 

Under the assumption of a simple radial flow, to explain 
the presence of multiple absorption line profiles that display 
complexities both in the velocity and ionization state parame- 
ter spaces, one typically is pressed to invoke a multiple-phase 
absorbing gas at different locations and with different veloc- 
ities. However, by relaxing the radial flow assumption and 
adopting a more realistic scenario, these complexities are nat- 
ural consequences of the geometrical extent of the outflow and 
proper treatment of cooling/heating processes. 

3.3. Effects of Unsteady Behavior 

So far, we have considered time-independent solutions 
of disk wind models. To illustrate the effects of unsteady 
solutions, we refer to the line-driven (LD) accretion disk 
wind model as presented by Proga & Kallmanj (120041) . These 
authors presented hydrodynamical simulations of a quasar, 
demonstrating how the accreting system is able to launch 
powerful equatorial disk winds. In particular, they found three 
different regimes for the flow: a 'hot flow' in the polar region, 
where the density and the velocity are very low; a 'fast stream' 
in the equatorial region, where a wind successfully developes, 
and a 'transitional zone' in between the two, where the flow 
is hot and struggles to escape the system. The 'transitional 
zone' is effectively shielding the 'fast stream' flow from the 
strong central ionizing continuum, preventing it from becom- 
ing over-ionized and losing acceleration in the UV resonant 
absorption lines. Broad-band X-ray spectra computed using 
detailed Monte Carlo simulatio ns to treat t he radi ative trans- 
fer in the flow are presented in I Sim et"an (1201 Obi) . Here we 
limit our study to the analysis of simple absorption profiles 
computed under the Sobolev approximation as in the previous 
Sections, using as input the flow physical parameters as given 
by the Proga & Kallman (2004) hydrodynamical simulations. 

In Figure [8] we show the absorption line profiles computed 
for five different inclination angles through the 'fast stream' 
and the 'transitional zone' of the flow, / - 60,65,70,75, and 
80° from top to bottom. Figure |9] shows the radial profiles of 
the velocity, density, and optical depth for the same inclina- 
tion angles as above. What is immediately clear is the depar- 
ture from the spherical symmetry: both the radial velocity and 
the density are not monotonic. Interesting physical insights 
can be obtained when looking at the 65° and 70° LOS, that 
are the red and green lines in Figure [8] and Figure |9] These 
two inclination angles roughly correspond to the transitional 
zone from the slow to the fast stream zone. Observationally, 
the 65° red spectral line profile shows much more blueshifted 
absorption than the 70° green profile. However, when looking 
at the radial velocity plot in Figure |9l it can be seen that the 
absorption at high velocity for / - 65° happens close to the 
origin (R ~ ISORq, where Rq - 8.8 x 10 cm is the inner 
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Figure 8. Blueshifted absorption line profiles computed for the LD accre- 
tion disk wind model of Proga & Kallman ( 2004) as seen by the different 
inclination angles i = 60, 65, 70, 75, and 80° from the disk rotational axis for 
the black, red, green, blue, and light blue, respectively. The velocity scale has 
been normaUzed to the maximum outflow velocity u„jax = 2.9 X lO' cm s"' . 

radius of the accretion disk and is equal to 3 Schwarzschild 
radii), where the dynamical instabilities of the transitional 
zone of the wind are the highest. At larger radii, for / - 65°, 
the radial velocity decreases: the deep absorption signature 
at Voui ~ 0.65i»,„flv is due to a portion of the wind (a 'puff', 
a blob) that will shortly (on dynamical time scale) fall back 
toward the disk plane. On the other hand, for / = 70°, the 
maximum of the absorption is reached at much lower veloci- 
ties, and shapes a broad deep trough at Vo„t ~ 0.1 - Q.3v,„ax- 
One can see that most of this absorption happens at large radii, 
R ~ 300/?() , and that this portion of the wind is more stable 
and effectively accelerated to leave the computational domain. 

Among the the considered LOS for this particular model, 
only the largest inclination angles (70, 75, and 80°) track por- 
tions of the flow that will successfully escape the system, al- 
beit they show the less blueshifted absorption profiles (Fig- 
ure [8]). On the contrary, the highly-blueshifted profile as- 
sociated to the 65° tracks a part of the flow where the gas 
is initially accelerated to the highest velocities, but being so 
close to the central continuum source it becomes quickly over- 
ionized, so loses driving force, and falls back to the disk plane. 
This lack of one-to-one correspondency between blueshift 
and terminal velocity of the wind is again contrary to the ex- 
pectations of a simple radial flow of constant velocity. Dy- 
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Figure 9. From top to bottom velocity, density, and optical depth profiles 
as a function of radius, for the same inclination angles color coded in Fig- 
ure \§\ The radius in in units of the inner disk radius, Rp = 3j?5 where 
Rs = IGMbhIc is the Schwarzschild radius (see lProga & Kallmanll2004l 
for details). 
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Figure 10. Map of the ratio between the mass outflow rate computed under 
the spherical symmetry assumption, and by using actual values from the nu- 
merical simulations of Proga & Kallman (2004). The axes are in units of the 
inner disk radius, Rq = ^Rs where Rs = IGMbhIc^ is the Schwarzschild 
radius. 



namical effects can be extremely important in shaping the ap- 
pearance of absorption Hne profiles, and should be taken into 
account when deriving physical properties such as mass and 
energy fluxes associated to the wind. 

All these considerations can be also visualized in FigurefTOl 
that shows the map of the ratio between the outflow mass flux 
computed under the spherical symmetric assumption and the 
actual outflow mass flux. The strong blueshifted absorption 
seen along the / = 65° LOS corresponds to the structure close 
to the origin, where the mass loss rate estimated under spher- 
ical symmetry is maximally overestimated (up to a factor of 
100!). For / - 70°, the discrepancy between the two mass loss 



rates is much smaller, up to an order of magnitude. These sim- 
ple examples illustrate how different (and sometimes counter- 
intuitively) is a realistic non-spherical, unsteady wind from a 
simple, homogeneous, and steady spherical case. 

4. DISCUSSION 

We have computed synthetic absorption lines against a con- 
tinuum point source for different simulations of accretion disk 
winds, with the goal of investigating the effects of the flow 
properties on the appearance of the line profiles. The used 
method is as simple and general as possible; scattered and 
emitted radiation is neglected. All the models considered 
her e are Compton-thin in th e studied LOS, except than for 
the lProga & KallrnanI ( 120041) line-driven disk wind. A proper 
treatment of Compton-thick LOS requires scattered radia- 
tion to be taken into account, as has been done by means of 
multidimensional Monte Carlo radiative transfer treatment by 
[slm et al. ( 2010b). Such computations are however extremely 
time-consuming, hence not easily applicable to the study of 
various geometries and LOS for different disk wind models. 

We found that the geometry and kinematic of the flow 
can imprint their effect on the blue edge of the absorption 
line (Section 3.1). In particular, the edge will appear sharp 
when the radial velocity distribution and the opacity distri- 
bution along the LOS are ~ constant within the portion of 
the flow that is in resonance with the continuum, i.e. where 
absorption occurs. The edge will instead appear smooth 
when the radial velocity and the opacity have steep profiles 
along the LOS. The extreme sharpness of the blue edge is 
in part due to the Sobolev approximation used in this work; 
in real situations, thermal and turbulent line broadening are 
expected to slightly smear the profiles. However, this does 
not affect the general conclusion of a different level of sharp- 
ness of the blue edge for different radial velocity and opac- 
ity distributions across the flow. We note how UV obser- 
vations of the rare iron low-ionization broad absorption line 
quasars (FeLoBAL QSOs, that show blueshifted absorption 
in the Fe ii and Fe iii transitions) usually reveal much sharper 
blue edges than for the more common high-ionization BAL 
QSOs cases (e.g.. Hall et al. 2002). This might imply a dif- 
ferent geometry and kinematics of the outflow for these two 
classes of quasars (see for exarn ple iLamy & HutsemekersI 
l2000HFaucher-Giguere et"ai]|20T2h . 

The ionization state of the gas has also important effects 
in shaping the absorption line profiles: non-radial solutions 
to the disk wind problem offered by hydrody namical simu- 
lations can be steady yet have a non-monotonic distribution 
of the ionization parameter along the flow and as such can 
produce complex shapes of the absorption lines (Section 3.2). 
In particular, due to the complex distributions of the den- 
sity and the radial velocity across the flow, a higher ioniza- 
tion state gas could be slower than a lower-ionization state 
gas. Other authors already discussed possible dependences 
between the velocity and the ionization state of the wind, 
in particular referring to t he observed X-ray warm a bsorber 
properties. For example, [Founds & VaughanI (1201 ll) report 
about a positive correlation between the ionization state and 
the radial velocity for the outflowing X-ray absorbing gas ob- 
served in NGC 4051, thus suggesting agreement with a radial 
outflow model. On the other hand, Rami rez et al.l (2005) sug- 
gest the existence of an anti-correlation between ^ and y,,,,, 
for the warm absorber observed against the X-ray continuum 
of NGC 3783: this observation, along with the asymmetry of 
the absorption Hne profiles, made the authors conclude that a 
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spherical outflow model can not reproduce the observations 
of NGC 3783, and rather two outflowing components mod- 
eled using a kinematical radiation-driven wind model are re- 
quired (see also Ramirez 2011). However, for the NGC 3783 
case, the authors assume that the resonance emission lines 
are narrow, and therefore do not affect the blue wing of the 
residual absorption features, while several examples of broad 
soft X-ray resonance em ission lines have been reported in 
flie literatoe (e.g ., Kaasti- a et al.ll2002l; ICostantini et al.ll2007l; 
iPounds & Vaughani201 li) . The situation is thus far from being 
clear, and only the launch of future X-ray missions capable of 
providing a high spectral resolution over a broad energy range 
(e.g., the microcalorimeter onboard ASTRO-H) will allow us 
to firmly assess the character of such correlations between the 
ionization state and the velocity of the X-ray absorbing gas. 

Not all known disk wind solutions are time independent. 
Therefore we also studied synthetic line profiles from highly 
dynamical models. We found the dynamics of the wind could 
strongly contribute to the observed complexity of the absorp- 
tion line profiles, especially along the LOS that track the most 
unsteady portions of the flow (Section 3.3). There may not be 
a direct proportionality between the observed velocity shift 
of the lines and the mass outflow rate, and the latter can de- 
pend in complex ways on the physics and on the geometry of 
the accretion disk wind. It appears that detailed simulations 
of synthetic spectra for accretion disk winds are necessary to 
provide realistic estimates of the mass outflow rate. One of the 
best studied examples is the X-ray absorbing wind observed 
in PG 1211-1-143 (Pounds et al. 2003). A mass outflow rate 
estimate using a spherically symmetric, constant velocity, ge- 
ometrically thin shell approximation for the interpretation of 
the high-ionization, high-veloci ty Fe xxv and Fe x xvi absorp- 
tion fines is Mom ~ 8.7Mq yr"' (iPounds & Reeves 2009), that 
after being corrected for a covering fraction of ~ 0.4 trans- 
lates to Moui ~ 3.5M0 yr"'. By fitting the same data to the 
models drawn from kinematic disk wind models and adopt- 
ing a proper treatment of the radiative transfer in the flow, the 
mass outflow rate estimate is found to be Mgui ~ 2.1 Mq yr"' 
jSim et al.ll2010al) . i.e. a factor ~ 40% lower. In the future it 
would be interesting to fit the data of large samples of AGN 
that display X-ray absorbing outflows to models de rived from 
hydrody nar nical simulations as those illustrated by lSim et al.l 
(1201 Obi) and lSim et aTl (120 12h . to derive refiable mass outflow 
rates. 

Distance estimates for the absorbers often rely on the as- 
sumption of a single-zone of gas photoionized by the central 
continuum source, and exploit the ionization parameter def- 
inition ^ = L/nR^ to get constraints on the location of the 
gas based on the observed density, luminosity, and ioniza- 
tion state. For low-density, low-ionization state gas this trans- 
lates to very large radial distances, of the order of the kpc 
(e.g., Korista et al. 2008; Moe et al. 2009; Dunn et al. 2010). 
Further assuming a constant velocity radial outflow, huge 
mass outflow rates are derived (up to several hundreds of so- 
lar masses per yea r). However, as already pointed out by 
lEverett et al.l (120021) . the derived radial location of the wind is 
strongly lowered in the case of a multi-phase absorber, where 
the low-ionization gas is shielded against the ionizing con- 
tinuum source. Multiple phases of absorbing gas are exactly 
what is observed. 

It is certainly convenient to think about AGN winds as 
spherically symmetric, constant density, geometrically thin 
shells of gas expanding with uniform velocity; we then as- 
sume photoionization equilibrium to derive constraints on the 



distance of the absorbing material from the continuum source. 
However, the location of the origin of the wind and of the 
origin of the absorption lines are two different quantities that 
may or may not coincide. For example, a wind can originate 
from a very small distance from the SMBH, yet an absorption 
line can be formed at larger radii, or over a large radial range 
(e.g. from a few tens up to several hundreds of gravitational 
radii). We have shown that if the wind does not resemble 
a 'thin shell' geometry but is rather radially extended, then 
the distribution of the kinematical, ionization, and dynamical 
properties of the flow along the LOS all come into play into 
shaping complex absorption profiles. Therefore, absorption 
lines from accretion disk winds form in regions that only par- 
tially cover the source: the covering fraction is dependent on 
the geometry, on the velocity, and on the ionization state of 
the gas (Proga et al. 2012). 

Complex and variable blueshifted absorption profiles are 
starting to be observed in the X-ray spectra of a num- 
ber of AGN; notably, PG 1126-041 and NGC 4051 both 
show variability of complex absorptio n structures on very 
short time scales in the iron K band dGiustini et al.l 1201 It 
IPounds & Vaughanll2012h . X-ray broad absorption lines are 
also observed to be strongly variab le in shape and strength 
(IChartas et al.ll2009t iLanzuisi et aLlf2012l) . AH these obser- 
vational results strongly point against a uniform and steady 
spherical outflow. 

5. CONCLUSIONS 

We have considered several physical models of accretion 
disk winds in AGN, and computed absorption line profiles 
predicted by such models. For the simplest isothermal disk 
winds, we found that although the models are non-spherical, 
they predict shapes of the synthetic profiles that could have 
very similar characteristics for different wind geometries. 
Therefore we conclude that the wind geometry is not the main 
contributor to the great diversity of the observed line shapes. 

More complexities and diversities of the profile shapes can 
be produced by facilitating one of the key properties of AGN 
disk winds: the non-monotonic distribution of basic wind 
properties such as velocity, density, and ionization state. Ow- 
ing to the complex distribution of these quantities along the 
LOS, multiple and detached absorption troughs can be easily 
produced in radially extended, non-spherical outflows. The 
highly dynamical nature of certain portions of AGN disk 
winds can also have significant effects on the mass and energy 
fluxes estimates, that can be off up to two orders of magnitude 
with respect to estimates based on a spherically symmetric, 
homogeneous, and constant velocity wind. 
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